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Microscopic cascading of second-order nonlinearities between two molecules has been proposed 
to yield an enhanced third-order molecular nonlinear-optical response. In this contribution, we in- 
vestigate the two-molecule cascaded second hyperpolarizability and show that it will never exceed 
the fundamental limit of a single molecule with the same number of electrons as the two-molecule 
system. We show the apparent divergence behavior of the cascading contribution to the second 
hyperpolarizability vanishes when properly taking into account the intermolecular interactions. Al- 
though cascading can never lead to a larger nonlinear-optical response than a single molecule, it 
provides alternative molecular design configurations for creating materials with large third-order 
susceptibilities that may be difficult to design into a single molecule. 



INTRODUCTION 



exceed the fundamental limit. 



Cascading is a process where photons generated by 
one molecule interact with nearby molecules. A lin- 
ear local field model is the simplest form of cascad- 
ing. When photons are generated through nonlinear 
processes that mediate interactions between molecules, 
more complex behavior is observed. The net effect is 
two molecules that interact through fields generated by 
a second-order nonlinear-optical process that act as a 
third-order nonlinear-optical process. Thus, it is natural 
to ask whether or not it is possible to make two-molecule 
cascaded systems with a larger nonlinear-optical response 
than that of a single molecule. [THi] 

There are two limiting geometries of aligned one- 
dimensional molecules that encompass all possible cas- 
cading susceptibilities. Previously, we discussed the cas- 
cading contribution of two side-by-side molecules with 
the long axis oriented in the direction of the electric 
field. [5] In this paper we discuss the effects of cascad- 
ing in a two-molecule system where the molecules are in 
an end-to-end configuration and aligned parallel to the 
applied electric field. 

Over the last decade, it has been shown that there 
are fundamental limits to the nonlinear molecular 
susceptibilities. [SHE] These limits have yet to be bro- 
ken. Moreover, experiments show a large gap between 
current measurements of the best second hyperpolariz- 
abilities and the theoretical limit. [TSHIS] Thus, it would 
be interesting if cascading could be used as a way to 
exceed the fundamental limits. Since the fundamental 
limits hold for any quantum system, a two-molecule cas- 
caded system should also obey the limits. However, since 
the intermediate photons can, in principle, be resonantly 
enhanced, it is not clear whether or not the limits can 
indeed be broken. This work seeks to better understand 
the nature of the cascading contribution to the second 
hyperpolarizability and to understand if it is possible to 



CLASSICAL ROUTE TO CASCADING 

Consider two freely-rotating one-dimensional 
molecules that are located on the same electric field line 
of a one-dimensional applied field, = EaZ^ where we 
are using the cartesian coordinates x, and z. When 
the molecules are also aligned with the field, the system 
is said to be in the end-to-end configuration. Since 
the molecules are assumed to be identical, |pi| = |pj|, 
where |p| represents the magnitude of the electric dipole 
moment. 

The induced dipole moment of the zth molecule in 
terms of the applied field and the electric field of the 
jth molecule is 

Pi = h ((/>i, Oi) a [/, ((/),, Si) • (E, + E,) 

+ h {^^. 0^) 7 [h {^^, Oi) • {Ea + E,)] ' • (1) 

Here, 

fi {4^1-, ^i) — sin(/)^ sin^^f — cos(/)^ sin^^^ + cos^^i, (2) 

where (j) is the azimuthal angle and 9 is the polar angle 
of molecule i. Note that i ^ j. Also, a, /3, and 7 are the 
polarizability, hyperpolarizability, and second hyperpo- 
larizability. The induced dipole field of the jth molecule 
at the iih. molecule's position in Equation [T] is 



E 



Pi 



J — 2 {—X sin sin Oj + y cos sin Oj + 2z cos Oj ) . 



(3) 

We wish to solve the self-consistent equation for 
Pi {6i^6j,(j)ij)^ where = — The contribution 
of the ith molecule to the first three effective molecular 



2 



susceptibilities, Ofeff,*, /^eff,*, and 7eff,i, are shown in the 
appendix, where they are derived by solving Equation [l] 
for Pi and using 



1 d^'p 



(4) 



where ^(^^ = a, ^^2) = and ^(^^ = 7. 

The polarization is maximum when the two molecules 
are aligned in the direction of the applied external electric 
field as shown in Figure [l] Thus, to study the upper 
bound of this case, we set Oi = 62 = 0, and 



Ei = E. 



2p 



(5) 



For the perfectly aligned case, the first three effective 
molecular susceptibilities, aeff, /^eff, and 7eff, yield 



aeff = 2r^ 



/3eff = 2r^ 



2a' 



/3 



(r3 - 2a) 



3 ' 



and 



7eff 2r 



12 



7 (r^ - 2a) + 4/32 
(r3 - 2a)^ 



(6) 
(7) 

(8) 



Note that the molecular susceptibilities. Equations [6]- 
[Sj are in scalar form. For typical 7r-conjugated donor- 
acceptor molecules, /3xxx = P dominates the response. 
All further calculations presented in this paper neglect 
the permanent molecular dipole moment and therefore 
the ground state transition moment is xqo = 0. Also, we 
will only consider the far off-resonant case, £^(cj = 0), 
where uo is the angular frequency. 
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FIG. 1. (Color online) The configuration of the molecules 
studied in this work. The two short black arrows represent the 
molecules while the long black arrows represent the applied 
electric field. The red (long dash) and blue (short dash) curves 
represent the electric field lines from each induced dipole. 



The polarizability has units of volume, and is loosely 
related to the 'size' of the molecule. The self-consistent 
calculation assumes two distinct molecules, so their sep- 
aration must be larger than their size. Therefore, the 
separation must obey r > \f2oL. Moreover, the factor, 

I (r^ — 2a), which can be viewed as a local field factor 
due to interaction, can be expanded to yield 



(9) 



2a 



The first term on the righthand side of Equation |9] de- 
scribes the classical cascading approximation with no lo- 
cal field, i.e. ^ 2a. In this approximation Equation [s] 
reduces to 



'eff 



27 + 8 



(10) 



where the second term on the righthand side is the cas- 
cading contribution. In the end-to-end geometry, the re- 
sultant cascading term is a factor of 8 larger than that 
given by Baev et a/. [4 

Before proceeding, we note that nonlinear optics is 
based on the series expansion. For convergence, each 
term must be smaller than the preceding one. When 
= 2a, the series does not converge, even for an in- 
finitesimal applied field, since /^eff/o^eff, 7eff//^eff, etc. di- 
verge. Thus, the fundamental basis for nonlinear optics 
fails when ^ 2+a. In this limit, the system is better 
treated as a single molecule, in which case the fundamen- 
tal limits must apply directly. 

The approach of this work is to test the regime of 
strong molecular interactions. If the fundamental lim- 
its are not exceeded in this regime, and are known not 
to be exceeded for two noninteracting molecules as well 
as for the single merged molecule, then this would prove 
that molecular interactions cannot lead to a second hy- 
perpolarizability that exceeds the limits. 

In the classically interacting approximation (r^ ^ 2a) , 
Equations |6][7| reduce to a polarizability and hyperpo- 
larizability of 2a and 2/3, as expected for two non- 
interacting molecules. However, there is an additional 
contribution to the second hyperpolarizability due to cas- 
cading as quantified by Equation [To) 



MOLECULAR INTERACTIONS 

The above classical model assumes that the molecules 
are unchanged when they are in close proximity to each 
other. However, we would expect interactions to cause 
the energy levels to shift and the transition dipole mo- 
ments to be altered. Since these effects get stronger as 
the molecular density increases, as do cascading effects, 
they need to be taken into account. 

We begin by considering the energy shifts due to molec- 
ular interactions. A system with n molecules can be ap- 
proximated by M"^ states, where M is the number of 



3 



states used to characterize a single molecule. For simplic- 
ity, we will approximate each molecule using a two- state 
model. The two-molecule system will then be described 
by four states. Thus, the resulting two-molecule system 
should yield a good approximation to the nonlinear re- 
sponse when it approaches the fundamental limit, [6 in 
which case E20 00, where Eij = Ei — Ej. While real 
systems are known to have an infinite number of states, 
numerical optimization studies of potential energy func- 
tions show that a quantum system, when at the limit, 
is dominated by three states. [Mj [16] Given that our in- 
vestigations focus on these limiting cases, the three-level 
model is appropriate. The fourth state provides a cor- 
rection factor when the hyperpolarizabilities are further 
from the limit. 

Two molecules in close proximity experience a van der 
Waals interaction, which in the end-to-end geometry is 
given by 



Vir) 



(11) 



where r is the distance between the two molecules, is 
the position operator for the zth molecule, and e is the 
charge of an electron in gaussian units. We use pertur- 
bation theory to describe this interaction. 

For a molecule with no permanent dipole moment, the 
perturbed energies for the geometric configuration shown 
in Figure [1] are 



^0 — ^0 
Ef = El 



Eq 
Eq 



■AE- 



1,0' 



where 



El 
El 



E1+E0 + AE+0, 
Ei + Ei + AEi,, 



'd^ E 



10 

g2 

^^r:o = -2^kioi , 

AP -9^^ 1^ |2 I kiol^ 



(12) 
(13) 
(14) 
(15) 



(16) 
(17) 
(18) 
(19) 



and the superscript 6* is a label for the coupled system. 
The transition moments Xij and energies Ej are for the 
unperturbed molecules. For a molecule with no perma- 
nent dipole moment, xqo = 0, xn can be replaced with 
Axio = ^11 — ^00 in the first term on the righthand side 
of Equation [19] - allowing the sum rules to be used to 



recast the above expressions. 

Next, we consider the effect of dipole interactions on 
the transition moments. The diagonalized wavefunctions 



for the four-level system are 
|0^) = |0a)|0b), 
|1^) = -^(|0^)|1b) 

|2^) = -^(|0^)|1b) 
|3^) = |U)|1b), 



|U>|Ob» 
|U>|Ob)) 



(20) 
(21) 

(22) 
(23) 



where the subset of degenerate states has been diagonal- 
ized with respect to the interaction Hamiltonian. Since 
there is no applied magnetic field, and thus no mag- 
netic vector potential in the conservative Hamiltonian, 
the wavefunctions can be chosen to yield real transi- 
tion moments. As a result, the transition moment ma- 
trix is symmetric, that is, Xij 
molecules are indistinguishable, x, 

The linear and nonlinear susceptibilities are calculated 
using a sum-over states expression. For the two-molecule 
system, the transition moments are calculated using the 
wavefunctions given by Equations |2Q]|23l the position op- 



Xji. Furthermore, the 

A ^ 

ij ^ij' 



erator x 



- X , and the energies given by Equations 
T2][T5 Using this model of the van der Waals interaction, 
the second hyperpolarizability is a simple sum of the form 
_ _^^B ^ direct second hyperpolarizability 

contributed by molecule i is given by 



n,m,l 
S\ -i \ jS\ 



X {m^\x' \l^){l^\x' |0^) 
(E^o) 

X {0^\x' \m^){m^\x' \0^) , 



(24) 



where the prime denotes that the sum excludes the 
ground state. It is straightforward to show that there 

Likewise, the hy- 



are no cross terms of the form x"^x^ 



TABLE L Transition moments 
Transition Moment Molecule A Molecule B 
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perpolarizability contributed by each molecule is 



I x" |m" ) (m" | x^ \{)^) 



, (25) 



where = + /3^. 

The barred position operator is defined as {i\x\j) = 
{i\ X \j) — (0| X |0) 5ij, where 5ij is the Kronecker delta. [T7] 
Thus, 

{i\x^ \j) = {kA\{lB\x^\mA)\nB) 

= {{kA\x^ \mA) - {Oa\x^ \0a) Skm) Sin- 

(26) 

The transition dipole moments for molecule A and B are 
shown in Table [H 

The transition energies used to calculate /3 and 7 from 



Equations 24 through 25 are derived from Equations 12 
T9I and are of the form, 



E 



10 



E. 



20 



^10 
^10 



^Foil +— ^Foil 
koil 



^30 — 2^10 



2e 



2e4 4 

-^FOll : 



I Ax 



4e4 



101 



Foil 



(27) 
(28) 
(29) 



CASCADING AND THE FUNDAMENTAL LIMIT 
Case 1: Quantum and classical interactions 

The generalized Thomas-Kuhn sum rules have been 
used to calculate the upper limit for the nonlinear molec- 
ular susceptibilities, [H [18] and are calculated without ap- 
proximation using the double commutation relationship 
between a conservative Hamiltonian and the position op- 
erator. The resulting sum rules are of the form. 



00 

E 

n=0 



1 



{Eg + Ep) Xq„X 



2m 



(30) 



and relate the matrix elements of the position operator 
to the energies. All solutions to the Schrodinger equation 
must obey the sum rules. 

The three-level ansatz states that when a quantum 
system has a nonlinear susceptibility that is near the 
fundamental limit, three states contribute. ^3 ^0] 
Thus, the four- level system resulting from two interacting 
molecules should be sufficient for studying the behavior 
of the largest cascaded nonlinear response. In the three- 
level ansatz, the normalized transition moment from the 
ground state to first excited state, X, is an important 
parameter, which is given by 



X 



Fiol 



(31) 



where xf^^ is determined from the sum rules to be of the 
form 



^10 



(32) 



2mEio 

Setting g = 0, we can multiply both sides of Equation 



30 by XpQ and sum over p, which gives 



|Xno| AXnO = 



00 ^ 



p^n 



nO 



^On^np ^pO' 



(33) 



Under the three- level ansatz, one can use the sum rules 
to find X12. Then, using Equation 33 and allowing £^20 ^ 
00, we get 



|xio| Axio 



-2Ko""lVl-X4. 



(34) 



Regardless of the opposite signs for some of the entries 
in Table |ll = and 7^ = 7^. To simplify the anal- 
ysis that follows, we define the dimensionless parameter. 



where 



mElo ' 



(35) 



(36) 



The individual molecular contributions, /3pert and 7pert7 
for the interacting system of two-level molecules with per- 
turbed energy states in terms of r' and X are calculated 



from Equation [25 

/^pert 



6\/2 



3/2^ 



7/2 
10 



(4r'i2 + X^y 



and Equation |24[ 



(37) 



Tpert 



nfiEl^ 



[(32r' 



/21 



4^/3^12^ (4 - gx"* + 

Sr'^X*^) (12 - 29X'' + ITX' 
■ -16) 
4) 



bX^ 



(4/12 ^ ^8^- ^2/^X2 + 4r'3 (1 - X^) + X 



(38) 



It should be noted that these solutions are based on 
perturbation calculations, which are only valid when the 
distance between molecules is sufficiently large to lead 
to small energy shifts. Since the electron cloud defines 
a physical size of a molecule and two molecules cannot 
overlap, a minimum possible separation is defined, as has 
been described for van der Waals forces between neigh- 
boring molecules. [2 H 12 2 j There are also retardation ef- 
fects which can lower the cascading contribution at large 
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FIG. 2. (Color online) The four- level quantum model of the 
second hyperpolarizability of two single molecules. When 
r' < \/2, the two molecules overlap and the calculations break 
down. 



distances. However, such effects are of no concern here 
since we are considering the maximum cascading contri- 
bution for two molecules in close proximity. 

The fundamental limit of the second hyperpolarizabil- 
ity is 



W\ 

^""f ^ 



(39) 



r3 
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FIG. 3. (Color online) A diagram illustrating the theoretical 
approach to cascading: (a) The relationship between 7^^^* and 
the individual contributions, (b) the second hyperpolarizabil- 
ity of two molecules each with a set number of electrons with 
respect to the second hyperpolarizability of a single molecule 
with the same number of electrons, and (c) a comparison be- 
tween the interacting maximum, 7max, and the noninteracting 

case, 7max- 



where the superscript 5* denotes the values of those pa- 
rameters for the two-molecule system. Note that the two- 
molecule system is characterized by a lower transition 
energy and has twice as many electrons as does a single 
molecule. The fundamental limit of the second hyper- 
polarizability for two molecules, each with N electrons, 
scales as {N^) = 4A/'^, thereby quadrupling the limit 
relative to the single molecule. 

The transition energy between the first excited state 
and the ground state becomes smaller as the two 
molecules are brought together, and therefore E^q < Eio 
- further increasing the limit. In the semi-classical ap- 
proximation, the basis of dielectric theory, the linear and 
nonlinear susceptibilities do not change in the presence of 
perturbations, thus, the eigenenergies of the underlying 
molecules are inviolate. The semi- classical fundamental 
limit, at r ^ 00 and X = 0, can then be expressed as a 
function of the number of electrons in each molecule. A/", 
and transition energy, £^10, of each molecule. 



7n 



16 



(40) 



Thus, the semi-classical approximation yields a limit that 
is smaller than the exact one, or 7max ^ 7max- There- 
fore, 7eff/7max ^ 7eff/7max and the semi-classical theory 
will give a larger intrinsic value than the exact result. 
If this ratio does not exceed unity in the semi-classical 
approximation, it will not do so for the exact case. 

Although 7max for specific values of X and r' is difficult 
to calculate analytically because Equation [38] is a high- 
degree polynomial of X and r', it can be numerically eval- 
uated. The numerical solution yields 7max ^ l-'^67max- 
However, this large value is found when r' < v^, the do- 
main where the separation between the molecules is less 
than their size - a physically unallowable regime. Figure 
2 1 shows a plot of 7sys = 7pert (^"^^^fo) using Equation 
33 for the two molecule system, where 7max and 7max are 
derived. Note that 7max — 7sys = 0, r ^ 00). If the 
intrinsic effective second hyperpolarizability in the semi- 
classical approximation, 7eff/7max5 of the coupled system 
never exceeds unity, then 7eff < 7max < 7max- If this 
is the case, then the effective second hyperpolarizability 
of two cascaded molecules will never exceed that of one 
molecule with twice the number of electrons. Figure [3] is 
a schematic diagram that illustrates the various terms in 
the cascading process and the the energy levels in each 
case. 

The effective second hyperpolarizability has a contri- 
bution from the hyperpolarizability due to the second 
term on the right hand side of Equation 10 Figure [4] 
)ws the two-level model for PpQ^t/^^^ term in Equation 
as a function of X and r' in units of e^h^N'^ /m?E\Q. 
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There is no cascading contribution at large distances, as 
expected. 

Substituting the perturbed system's values for 7pert 
and /3pert from Equations 37 and 38 into Equation 10 



6 



gives the effective second hyperpolarizability with cor- 
rected transition energies. Figure [5] is a surface plot of 
the effective second hyperpolarizabihty with perturbed 
energies divided by 7max- As r ^ oo the intrinsic hy- 
perpolarizabihty converges to 1/2, as is expected for the 
limit of a two electron molecule, 7sys, having {N^) — 
very large distances, the two molecules 
are independent with 7 = 7^ 4-7^ oc 2N^ . Note that the 
red plateau region in Figure [5] shows the domain where 
the limit would have been exceeded. 

The second hyperpolarizability never exceeds the fun- 
damental limit for ^ 2+. This is an important consid- 
eration for molecular design because increasing the po- 
larizability increases the minimum required separation 
distance, thereby reducing the cascading contribution. 
Moreover, the limit is never surpassed for any distance 
that is greater than \l2oL. This length is of critical impor- 
tant because any molecular separation that is less than 
this distance is unphysical and can lead to a divergence 
of the local field factor. 



Case 2: Only classical interactions 

While a quantum characterization of interactions be- 
tween two molecules leads to interesting insights, it is too 
complex to solve. In this section, we consider the special 
case where the energies are unchanged by the interaction, 
and we use the classical cascading approximation given 



by Equation 10 



Due to the simplicity of the above approximations, we 
choose to describe each molecule using three states. We 




3 1.0 



FIG. 4. (Color online) The cascading contribution of 
Ppevt/^'^ • The perturbed Hamiltonian keeps the function from 
diverging as r' 0. 



define the fractional transition energy, as 

^10 



E 



E' 



20 



(41) 



When E20 = Eio^ E = 1^ and when E20 cx), E = 0. 
Thus, the system's energy spacing is fully parameterized 
by E in the range from zero to unity. 

The sum rules for a three level model have been shown 
to give [6^ 



1^121 



^10 



1-E 



(42) 



Substituting Equation [42] into Equation [33| gives 

yi-x4, (43) 



l^iol Axio 



,2^-2 



^10 



which is the three level generalization of Equation [34] 
This can be understood by letting £^20 ^ 00, thereby 



reducing Equation [43] to Equation [34] Similarly, for the 

second state, we get 



|x2o| Ax2o = Ko^l' (1 - 2E) ^ ^A_xvT+X^. 

(44) 

The hyperpolarizability and the second hyperpolariz- 
ability of a single molecule with no external infiuences 
can then be written as 



I3 = 3e 



3 / ^On-^nrri'^mO 



(45) 



y(0) 




FIG. 5. (Color online) The intrinsic effective second hyperpo- 
larizability of two interacting molecules. The thick black line 
is the lower limit of the center of mass displacement given by 
the condition r' > \/2. The thick red line is the experimental 
limit of closest approach for spherical atoms and molecules 
with large polarizabilities, r > 2. [23] The values that exceed 
the semi-classical limit have been removed from the diagram, 
and are shown as the red plateau. 
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FIG. 6. (Color online) The top two graphs are the intrinsic 
hyperpolarizability and second hyperpolarizability as a func- 
tion of X and E for a single molecule system. The bottom 
four graphs show the intrinsic effective second hyperpolariz- 
ability of the aligned molecules with no energy corrections as a 
function of X and E for several r' values. The {3 contribution 
decreases as the distance of separation increases. 



and 



n,m,l 
oo 



(46) 



These sum-over-states expressions were first derived us- 
ing quantum perturbation methods by Orr and Ward. [17] 

Substituting Equations [3lJ |32) and [iTpi] into Equa- 
tions [45] and [46] we find the unperturbed three level ex- 
pressions for P and 7 in terms of X and which are 



3 e^h^ 

l^-E^E'^]xVl-X^ 



^3L 



m J 



(47) 



'eff 
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FIG. 7. (Color online) The semi-classical intrinsic effective 
second hyperpolarizability of two end-to-end molecules. The 
thick black line is the minimum separation distance given by 
the electromagnetic size of the molecule, or r' > \/^. The 
thick red line is the experimental limit of closest approach for 
spherical atoms which gives r > 2.(23] The red plateau is the 



region where 7eff /7n 



< 1. 



and 



73L 



m?ElQ 



[4- 



lf{E 



2{E^ - 1)E^X^ 
1){E^ 



E + DX'' 



E + 3) E^] 



(48) 



The classically interacting effective second hyperpolar- 



izability is found by substituting Equations [47] and [48 
into Equation [lO] Figure [6] shows the intrinsic effec- 
tive hyperpolarizability for various separation distances. 
The second hyperpolarizability at the minimum separa- 
tion distance, r' = v^, never exceeds the fundamental 
limit. For a single molecule, the intrinsic second hyper- 
polarizability is bounded by —1/4 < 7^"* < 1. 

To illustrate how the effective second hyperpolarizabil- 
ity behaves as a function of r', we set £^ = so that each 
molecule is described by a two-level model that obeys the 
sum rules. The off-resonant nonlinear molecular suscep- 
tibilities are maximized in the two-level limit as shown 
in the top two graphs of Figure [6] This case models the 
effective second hyperpolarizability as a function of the 
separation distance for two interacting molecules with 
large nonlinearities. 

In the absence of a second excited state, the expres- 



sions given by Equations 47 and 48 reduce to 



^2L 



3/2 



X^J\ - X^ 



(49) 
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and 



72L 



(4 - 



(50) 



respectively, where the 2L subscript represents the sys- 
tem for £^20 -> oo. In this case of classically interacting 
molecules, there are no shifts in the energies of the eigen- 
states. For unperturbed molecules, 7 ^ oo for ^ - 
clearly an unphysical result when the two molecules are 
separated by a distance that is less than there electro- 
magnetic diameter. Therefore, we will restrict our dis- 
cussion to the intrinsic value of molecules with the same 
number of electrons as two separate molecules and re- 
strict the space to the physically-allowed region defined 
by r' > v^. 

Figure [7| shows the semi-classical intrinsic effective 
second hyperpolarizability, 7eff/7max5 which diverges as 
r ^ 0. This is a consequence of the classical model which 
assumes that the eigenenergies do not change even when 
the molecules interact with each other. Although the 
function diverges at zero separation, the second hyper- 
polarizability due to cascading never exceeds the funda- 
mental limit when the separation distance is larger than 
the electromagnetic size, r > \/2amax- 



REAL SYSTEMS 

In this paper and in the companion article, we have 
considered cascading only between molecules with fixed 
orientation. The degrees of freedom offered by cascading 
include the separation between molecules and their ori- 
entations. The average distance between molecules, the 
parameter studied by Baev et al^ can be easily controlled 
by adjusting the density of molecules in solution. The ad- 
ditional degree of freedom originating in the orientations 
of the molecules can be controlled by the application of 
an external electric field, which will also affect the lo- 
cal pair orientational correlation function due to induced 
dipole moments. Alternatively, local molecular forces, as 
are found in liquid crystals, could also have alignment- 
mediated enhancements. These real systems will lead to 
a cascading contribution that is lower than for an ideal 
fixed configuration. 

To understand how an externally applied electric field 
can be used to control cascading, we develop a thermo- 



dynamic model of poling that includes the effects of the 
local reaction fields as treated with the self-consistent 
field model. In a gas, the molecules will relax when the 
external field is turned off, but in a poled polymer, |24] 
the alignment can be made permanent. 

An external electric field applied to an ensemble of 
molecules in thermal equilibrium will result in a distri- 
bution of orientations centered along the field direction. 
External stress can be used to further control the orienta- 
tions, and therefore, the tensor properties of the nonlin- 
ear optical response. ^3 Here, we will consider the effects 
of an applied electric field only, and imagine that the ori- 
entational order is frozen in place as it would be in the 
case of a dye-doped polymer. 

We can understand the physics of the problem as fol- 
lows. The external field tends to align the molecules. 
However, the induced dipole moment of each molecule 
adds to the field at the site of the other molecule, 
where the local field at each molecule must be solved 
self-consistently. Thus, the orientational order of each 
molecule will depend on both the applied electric field as 
well as the distance between the molecules. We label this 
as the correlated case. 

If interactions between the induced dipole moments 
of each pair of molecules are not taken into account, as 
is the case in early models of poling, cascading effects 
will be enhanced, but to a smaller degree. The difference 
between the two cases is a measure of the local correlation 
strength. In each case, we assume that the orientational 
order has been frozen in place by the poling field after 
it has been turned off, and that the subsequently-applied 
optical field, which we approximate in the zero frequency 
limit to describe the off-resonance response, probes this 
orientational distribution. 

It is worth mentioning the differences between our ap- 
proach and that of Baev et al We are treating the end- 
to-end case, which our self consistent field calculations 
show to yield the largest response. In contrast, Baev et 
al use the side-by-side geometry. Additionally, we take 
into account the distribution of orientations using a ther- 
modynamic model rather than using a static orientation. 

Consider a pair of para-Nitroaniline (pNA) molecules 
in an ensemble, which are free to rotate in three di- 
mensions. The angular probability density for the ith 
molecule, P^, at an azimuthal angle, and polar angle, 
Oi, from Boltzmann statistics is 



Pi {(t)i,6i,(l)j,6j) 



lo C lo C d^^dOd^dO, sin 0, sin O.e^o Mdw^^d^)/kT ' 
I 



(51) 



where k is the Boltzmann constant, T is the temper- ature, Ej is the electric field due to the jth molecule. 
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6 is the polar angle measured from the z-axis, and (j) 
is the azimuthal angle. Note that for two molecules, 
1 and 2, there are only two possible combinations, i.e. 



(i,j) = (l,2) and (i,j) = (2,l). 

The exponents in Equation |5l] can be expanded in a 
power series to fourth order in the applied electric field by 
using the molecule's effective molecular susceptibilities, 



f 

Jo 



Pi • (d-Ea^d-Ej) 



L 



llEa COS Oi + ^aeff,i {(t)ij,0i,0j)El COS Oi + ^/3eff,i 



^- , 6>i , Oj ) El cos Oi + -7eff,i {(t)ij , Oi , Oj ) El cos Oi , 



(52) 



52 



orinate from • fi 



where the cosines in Equation 

Ea COS Oi, p^^ is the effective induced dipole moment that 
takes into account interactions between molecules using 
self-constant fields, and the permanent dipole moment, 
/i = 6.2 X 10~^^ erg^/^cm^ [26 is now included to accu- 
rately describe the correlation fields for pNA. To find the 
contribution of microscopic cascading, the denominators 



r 



in the zth molecule's effective susceptibilities, which act 
as the local field factors due to the polarizability of the 
other molecule, are expanded in a series and only the first 
term is kept. 

The orient at ion- averaged contribution to the zth 
molecule's effective second hyperpolarizability, (7eff,i), is 
then given by 



(7< 



JO 



27r 



27r 



d(j)idOid(j)jdOj sin Oi sin Oj cos OiPi 



/)7eff,i I 



(53) 



Furthermore, (7eff,i) = (7eff,2) due to the fact 
that the molecules are indistinguishable. Thus, the 
orientationally-averaged effective second hyperpolariz- 
ability due to the two- molecule system is (7eff) = 

2 (7eff,i). 

Para-Nitroaniline (pNA) is a donor- acceptor molecule 
that has been studied in great detail, and it is often used 
to test theoretical models. [26H28] It can be approximated 
as a one-dimensional molecule due to its strong push-pull 
behavior from the opposing ends of the benzene group. 
This molecule has previously been the subject of density 
functional theory calculations, where the second hyper- 
polarizability was reported to be enhanced due to micro- 
scopic cascading in the side-by-side geometry. [4] Here, we 
test the case of the effective second hyperpolarizability in 
the end-to-end geometric configuration in the presence of 
the poling field. 

The three-level ansatz will give accurate results for 
molecules with susceptibilities near the limit. [20 Since 
pNA is far from the limit {f3 and 7 are each about 2 orders 
of magnitude below the fundamental limit), we will use 
the measured values of the linear and nonlinear molecular 
susceptibilities instead of the calculated ones; but, will ig- 
nore the changes in the transition energies that arise from 
molecular interactions. Figures |5] and [7| show the regions 
where the classical model of cascading is a fair approxi- 
mation, that is, r' > \/2. The susceptibilities we use here 
are, a = 1.7 x IQ-^^ cm^ /3 = 9.2 x 10"^^ erg-^/^cm"^. 



and 7 = 1.5 X 10~^^ erg~^cm^ - values that were reported 
by Cheng et a/. [26 

Equations [57]|59| in the appendix show the effective po- 
larizability, hyperpolarizability and second hyperpolariz- 
ability of molecule i in the presence of molecule j for two 
arbitrary but fixed orientations, solved self-consistently 
in analogy to the simple end-to-end case given by the 
derivation leading to Equations |6] through [8| The orien- 
tational averages of these results are calculated accord- 
ing to Equation [53] These calculations, which are far too 
complex to solve analytically, were evaluated numerically 
using Mat lab.® 

First, we calculate the cascading contributions for the 
uncorrelated angular distribution for the tensor compo- 
nents along the applied field. For a molecule that makes 
an angle with the field, the effective polarizability and 
(second) hyperpolarizabilities along the field are given by. 



ouncorr 



= acos^, 
/3cos^ i9. 



and 



7effr =7cos3 



+ 4 



/32 cos4 e 



(54) 
(55) 

(56) 



The subscript i has been dropped from the angle because 
both molecules are identical and have the same orienta- 
tional distribution function due to the applied electric 



field. Then, the orientational average of Equation 56 is 
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FIG. 8. The effective second hyperpolarizability of two inter- 
acting pNA molecules at a temperature of 297 K where the 
separation vector between the two molecules is along the di- 
rection of the applied field. The molecules are freely rotating, 
where the alignment is a result of the applied electric field in 
units of StatV/cm The solid curves include interactions be- 
tween polarized molecules and the dotted curves include only 
the effects of the applied electric field. 



FIG. 9. The contribution to the effective second hyperpolar- 
izability of one freely rotating molecule that is coupled to a 
second molecule. The second molecule is held in both the di- 
rection perpendicular and parallel to the applied field while 
the other is free to move. The temperature is 297 K and the 
poling field is 10^ StatV/cm. 



determined using the orientational distribution function 
that includes only the effects of the applied poling elec- 
tric field. The correlated case is too complex to express, 
and is calculated using the self-consistent field calcula- 
tion to determine the orientational distribution function 
from which the orientational average of the effective sec- 
ond hyperpolarizability is calculated. 

Figure [8] shows the orientational- averaged effective sec- 
ond hyperpolarizability, (7eff), for the correlated (solid 
curves) and uncorrelated cases (dotted curves) as a func- 
tion of the distance of separation, r, for several electric 
field strengths. At large separation and zero applied field, 
the orientational average of the effective second hyperpo- 
larizability is (7e//) = "^"ixxxxl^i as expected. 

The applied poling field at a field strength of Ea = 10^ 
at 6A separation (the approximate size of the molecule) 
enhances the effective response, which includes cascad- 
ing, by appoximately a factor of four while correlations 
due to induced dipole interactions leads to another factor 
of about two. Thus poling correlations add constructively 
to the orientational ordering effects of the field, leading 
to approximately an order-of magnitude enhancement. 
Note this this is far short of the largest possible enhance- 
ments for perfectly aligned molecules due to thermal dis- 
ordering effects. 

In this example, while the net enhancement effect of 
the intrinsic hyperpolarizability is about an order of mag- 
nitude relative to the sum of the nonlinearities of two 
individual molecules, the largest possible nonlinear re- 
sponse when the two molecules are combined together 



into a single molecule would give a gain of a factor of 
2. Thus, the enhancement effect due to cascading yields 
an increase in the intrinsic hyperpolarizability by about 
a factor of 5. Thus, while cascading may not break the 
limit of the nonlinear-optical response, it can lead to an 
increase in the nonlinear response over the best avail- 
able molecules for applied fields that are at the dielectric 
breakdown limits for the best materials. 

The coupling strength between two molecules depends 
on the distance between them, and is a key parameter for 
cascading. The effect of the strength of interaction on the 
orientation can be isolated by allowing one molecule to 
rotate freely and the other to be fixed with its long axis 
either perpendicular or parallel to the applied field. Fig- 
ure [9] shows (7eff,i) for these two orientations of the fixed 
molecule. The inset, which shows a magnified view of the 
effective second hyperpolarizability of the perpendicular 
orientation, illustrates the weakness of the intermolecular 
interaction from orientational averaging in the presence 
of a strong applied fields when the fixed molecule is held 
perpendicular to it. 

These results give insights into cascading and how ex- 
tra degrees of freedom can be used to design molecu- 
lar systems with a larger intrinsic nonlinear-optical re- 
sponse. It is clear that cascading of lower-order nonlinear 
molecules with additional control of orientational order, 
as would be the case for liquid crystals or mechanically- 
stressed polymers, would further increase the microscopic 
cascading contribution due to the control of additional 
order parameters. 
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APPROXIMATIONS 

In this section, we briefly discuss the approximations 
used in our calculations, and how they may limit the 
validity of the results. 

Off-resonance response. All of the calculations con- 
sider the zero-frequency limit of the electronic response. 
In real experiments, the optical field frequencies can be 
tuned to be smaller than transition frequencies within 
the molecule, but there will always be some resonance 
enhancement. These can be partially taken into account 
by using the two level dispersion model by Gudar and 
Chemla.[29] 

Dipole approximation: Our calculations assume 
that the induced dipole moment can be expanded in a 
series of the electric field of the form p = aE + PE'^ + .... 
Many nonlinear-optical processes originate under the as- 
sumption that the series converges. This implicitly as- 
sumes that the fields must be small, that is, ^E/ j3 <C 1, 
etc. to all orders. Thus, the molecules must not get close 
enough together to violate this condition. 

Finite physical size: It is assumed that the sepa- 
ration between the molecules must be larger than their 
size, otherwise the two molecules will merge together to 
form a new one. The electromagnetic size of a molecule 
is determined from the polarizability a to be r > (2a)^/^. 

Three- level ansatz: The three- level ansatz states 
that when a molecule has a nonlinear response near the 
fundamental limit, it can be modeled as a three-state sys- 
tem using a sum-rule-restricted reduced-parameter space. 
While this has not been proven rigourously, it is an empir- 
ical observation that holds for a wide rage of systems. [20 
Thus, results in the regime of small nonlinear response 
are inaccurate; but, near the limit - the focus of this 
paper, the calculations are accurate. 

Magnetic interactions: We assume that the electro- 
magnetic vector potential can be ignored, so that that 
the position matrix is real. As a consequence, Xij = Xji. 
This assumption does not change the generality of the 
results. [30 

Perturbation approximation: The energy shifts 
and change in dipole matrix elements are calculated us- 
ing perturbation theory; so, the molecules must be far 
enough apart to limit interaction energies to be smaller 
than the eigenenergies of each molecule. Details of why 
this approximation holds at all distances is described in 
the companion article. [5 

pNA In the calculation of the self-consistent field 
interactions between two pNA molecules, the energy 
shifts and dipole matrix changes are neglected. This ap- 
proximation will not hold when the separation between 
molecules gets small, so the calculated results will not be 
accurate in the limit of small separations. 

ID molecules: This approximation is commonly used 
for push-pull molecules and does not lead to large devi- 
ations from the non-idealized case. The uncertainty in- 



troduced by this approximation for a system with two 
nonzero tensor components is about ^x/ + For 

example, if the next larger tensor component is 1/3 of 
the largest one, this will introduce an error of 5%. 

See the companion article for a more detailed discus- 
sion of the physics behind these approximations. [5] 



CONCLUSION 

In this work, we have investigated cascading as a 
method for enhancing the second hyperpolarizability. 
The fundamental limit has been used to define the scale- 
invariant intrinsic hyperpolarizabilities, which enables us 
to assess the effectiveness of cascading as a means of pro- 
ducing nonlinear susceptibilities larger than what is pos- 
sible for a single molecule. 

In the case of two one-dimensional molecules in an end- 
to-end configuration aligned with the direction of the ap- 
plied electric field, the cascading contribution to the sec- 
ond hyperpolarizability was shown to be bounded by the 
fundamental limit provided that > 2amax, a require- 
ment that the two molecules do not overlap. Addition- 
ally, this condition ensures that the fundamental basis 
of nonlinear optics, i.e. that the polarization can be ex- 
pressed as a power series of the electric field, is valid by 
ensuring convergence. 

Our full theory uses the self-consistent field approach 
to derive the effective nonlinear susceptibilities of the 
two-molecule system, and includes the effects of van der 
Waals interactions on the energies and transition mo- 
ments of each molecule. This full quantum treatment 
shows that the two-molecule system is bounded by the 
same fundamental limit as one molecule, the same re- 
sult as we get in the other extreme case of two side-by- 
side molecules. [5 Thus, we conclude that combining two 
molecules into one or allowing two molecules to interact 
with cascading leads to the same optimized response. 

The bottom line is that cascading is not a loophole 
that allows the fundamental limits to be exceeded. How- 
ever, cascading provides an additional degree of freedom 
for designing materials, and therefore may provide an al- 
ternative approach for improving the nonlinear response 
over what is possible through synthetic means alone. As 
we show for the model system pNA, cascading can yield 
an enhancement of the intrinsic second hyperpolarizabil- 
ity, i.e. the response can get closer to the limit. It may 
be difficult to manipulate nonlinear responses with a sin- 
gle molecule design of the same size as the two-molecule 
system that interacts with cascading. 
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APPENDIX: EFFECTIVE MOLECULAR SUSCEPTIBILITIES FOR TWO FREELY-ROTATING BODIES 



For two molecules of arbitrary fixed orientations and separation vector along the applied electric field, the effective 
polarizability of molecule i is derived from Equations fTE) and yields 



Q^eff,i = r a- 



(r^ + a cos {26 j) + a) cos 6i — a cos (pij sin Oi sin Oj cos Oj 
o? [cos (^ij (sin (20 i) sin (2Qj) — cos (^ij sin^ sin^ ^j) — 4 cos^ cos^ ^j] 



(57) 



where the local field factor is automatically and naturally taken into account through the self-consistent field calcu- 
lation. The effective polarizability of the two molecules together is aeff = Q^eff,i + <^eff,2- The next term given by 
Equation [4] is the zth molecule's contribution to the effective hyperpolarizability, which gives 



= [4r^a cos^ cos di - Ar^a cos Oj cos^ Oi + 2a cos^ Oj (4r^ cos^ Oi + 4a cos"^ 6>i - cos ^i^- sin Oj sin 6>i 
-f 2 cos^ Oi (r^ + a cos (j)ij sin (—2a cos Oi sin (26>j ) + sin Oj (r^ + a cos sin 6>j sin Oi))) 
— r^a cos (j)ij sin (26>j) sin (26>^)] 

/ 2 (r^ + 2a cos 6>j cos Oi — a cos sin Oj sin ^ (r^ — 2a cos Oj cos + a cos sin sin Oi) 



(58) 



The effective second hyperpolarizability is found in the same way. The ith molecule's contribution to the effective 
second hyperpolarizability of a system of two fixed molecules in the end-to-end configuration is given by, 

7eff,i = I r^a/3^ (—2 cos Oi cos Oj + cos (f)ij sin Oi sin Oj)'^ (cos Oj (r^ + a cos (2^^) + a) — a cos {Oi) cos sin sin Oj) 
X (4r^a cos Oi cos^ 6>j — 4r^a cos^ Oi cos 6>j — 4a^ cos (j)ij cos^ 6^^ sin Oi sin (26>j) 

r.2 /3 . /"/I^S ^^^2 /3 I /I ^, ^^,-,4 /3 '^'^ 'JUS (/^^j Sill (7^ Sill {7j j 

'^^asin {20 i) sin (2^^) cos ^ij)] 



-f 2a cos^ Oj (4r^ cos^ Oi + 4a cos^ Oi — cos sin Oi sin 6>j 
-f 2 cos^ (r^ + r^a cos (j)ij sin sin Oj + a^ cos^ 

/ (r^ + 2a cos Oj cos Oi — a cos 0^ j sin sin 6>j ) ^ (^^ — 2a cos Oi cos + a cos j sin Oi sin 6>j ) 



sin^ Oi sin^ ( 



r^7 (cos Oi (r^ + a cos {20 j) + a) — a cos 



sm Ui sm t/o- cos t 



(r^ — 4a^ cos^ Oi cos^ — a^ cos^ (j)ij sin^ sin^ Oj + a^ cos sin {20i) sin (26>j)) 
^a7 (2 cos Oi cos 6>j — cos (j)ij sin sin Oj) (cos (r^ + a cos (2^^) + a) — a cos (j)ij cos sin sin Oj) 



(r^ — 4a^ cos^ Oi cos^ 6>j — a^ cos^ sin^ sin^ Oj 



i(j)ij sin (2^i) sin (2^^))" 



+ (2 cos Oi cos 6>j — cos (j)ij sin 6^^ sin Oj) (2 cos (r^ + a cos (26>j) + a) — a cos (j)ij sin (26>j) sin Oi 

X (l6a^ cos"^ 6>j cos^ {Oi) + 8r^a cos^ Oi cos 6>j + 4r^ cos^ Oj (r^ + 4a cos^ Oi ^ a cos sin Oi sin 6>j) 

— 8a cos^ Oj (r^ cos 6^^ + a cos (j)ij sin (26^^) sin Oj) 

— a cos (4r^ cos^ 6>i sin sin Oj + sin (26>j) (2r^ sin (26>^) — a cos sin (26>j) sin^ ^i)))] 

I 4 (r^ + 2a cos Oi cos 6>j — a cos 0^ j sin sin Oj )^ [f* — 2a cos cos Oj + a cos (j)ij sin sin 6>j ) ^ j | 



/ 



1 



(—2 cos Oi cos + cos (pij sin sin ) 



(59) 



These results are used extensively is the section on real molecules. Equations [57| through [59) which describe the 
effective molecular susceptibilities of two interacting molecules, are exact for a uniform applied electric field and a 
nonlinear-optical response that can be represented as a power series in the electric field vector. 
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